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Abstract. The thermoelectric coefficient 1]{T) in unconventional superconductors is enhanced below Tc by intermediate 
strength impurity scattering that is intrinsically particle-hole asymmetric. We compute T]{T) for a strong-coupling d-wave 
superconductor and investigate the effects of inelastic scattering originating from electron-boson interactions. We show that 
Tj(r) is severely suppressed at temperatures just below Tc by a particle-hole symmetric inelastic scattering rate. At lower 
temperatures inelastic scattering is frozen out and T]{T) recovers and regains its large amplitude. In the limit T ^ 0, we have 
Tj(r) ~ r]^)T + ff[r^\, where the slope T70 contains information about the Drude plasma frequency, the details of impurity 
scattering, and the change in effective mass by electron-boson interactions. In this limit 1]{T) can be used as a probe, 
complementary to the universal heat and charge conductivities, in investigations of the nature of nodal quasiparticles. 



Low-temperature transport measurements have pro- 
vided a wealth of information about nodal quasiparticles 
in unconventional superconductors Thermal con- 

ductivity is of particular importance, because theory pre- 
dicts universality in the sense that the low-temperature 
asymptotic does not depend on the properties of the im- 
purity potential |3]. This prediction has been confirmed 
experimentally 1 1 ] . However, there are some difficulties 
in analysing the low-temperature thermal conductivity. 
First, the leading T^-dependence of the electronic con- 
tribution to k{T)/T is masked by a phonon contribution 
with the same T^-dependence |2]. Second, experiments 
are often done on ultra-clean samples, which means that 
the T^ power law of k{T)/T holds only in a small tem- 
perature bracket 14). 

In a recent paper we discussed how isotropic elas- 
tic scattering by impurities of intermediate strength, i.e 
described by a phase shift < < f , gives rise to 
an electron-hole asymmetric scattering time and con- 
sequently a large non-universal thermoelectric response 
01 ■ A careful study of the thermoelectric coefficient 
in the low-temperature regime would reveal informa- 
tion about the bare elastic scattering rate and potential 
strength. This information is hard to extract from thermal 
conductivity data. In this report we examine the interplay 
between elastic and inelastic scattering and show how it 
affects the thermoelectric coefficient. 

The thermoelectric coefficient, T7(r), is defined as 
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The quasiclassical propagator 5g^ has a closed form in 
which the self-consistently computed equilibrium Green 
function and the self-energy E^ serve as input (see 
Graf et al. 0] and Ref. |5] for details). In the present 
study we consider a composite self-energy in particle- 
E^, where the impurity self- 

and E«. 

its particle-hole symmetric and anti-symmetric parts, re- 
spectively. The self-energy E^ includes the effects of in- 
elastic electron-boson scattering, but we also assume that 
this interaction mediates the pairing. Below, is the 
off-diagonal component, i.e the usual strong-coupling 
function related to the energy dependent gap as A*(£) = 
W^{e)/Z{e). Contrary to the impurity self-energy the di- 
agonal part of the inelastic self-energy is particle-hole 
symmetric as Eq^^, = 0. Finally, Z(e) is the energy- 
renormalization function defined by the scattering renor- 
maUzed energy e'^ = Z*(£) £ = e - Ef ; - Ef This 
model is the same as we used in our study of the ther- 
mal conductivity fT). With this input we use Ref. jj] and 
write down the response function 



dee sech 



4r2j 2T 

^(p/,e)3E« (e 



dPf[vf,iVf,.i\ 



(2) 



where g) = —o/(^3 [e*/i2*(p/;e)] is the density 

of states, and = ^|W^(p/,e)|2 - (e«)2. Note that 
7] (r) is directly proportional to the imaginary part of the 
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FIGURE 1. Thermoelectric coefficient, Ti{T), calculated for four orders of magnitude in scattering rate, as indicated in the right 
panel. The dashed lines are results of a self-consistent calculation with the inelastic scattering model introduced in Ref. |4]. The 
full lines are the corresponding results obtained with an effective model where the inelastic part is included only as an effective 
mass via m* /m = 1 + Aj^ [4]. The temperature dependence obtained with the effective model is the same as within weak-coupling 
theory. In the left panel Tj {T) is scaled by the normal state Drude conductance, oxi = e^v^^Af-/2r, where vf, jVf, and T are effective 
values of respective bare constants. The arrows indicate the cross-over, T* , between the two scattering modes. In the right panel we 
show the temperature dependence of r]{T)/T with emphasis on the low-temperature regime where the leading T— dependence is 
r\{T)jT ^r\Q. The constant, t)o, is non-universal since it explicitly depends on the ratio of c?3Zq ,jjj^^(e) /5e|e=o and SEf ,„,^(e = 0). 



particle-hole asymmetric part Z^Qi,„p, which is an odd 
function of energy. 

In Figure n we show rj(r) for a strong-coupling d- 
wave superconductor. The results are compared with an 
effective model-calculation where only the mass renor- 
malization by inelastic scattering is accounted for, as de- 
scribed in Ref. 1^. As seen in Figure ^ inelastic scat- 
tering affects TiiT) much in the same way as it does the 
thermal conductivity reported in Ref. |4], /.e it dominates 
at temperatures T < but freezes out at low T <C Tc 
where instead elastic scattering limits the transport. We 
separate the two regimes where the two types of scat- 
tering dominate by introducing a cross-over temperature 
T*, defined as the temperature where ri{T) has its max- 
imium. The value of T* depends on the elastic scattering 
rate F for a given electron-boson coupling spectra. The 

suppression of rj(r) in the interval T*^T < Tc comes 
from that in our model inelastic scattering does not break 
particle-symmetry and hence reduces the weight of the 
kernel in the integral in Eq. (|2ji. 

In the low-T regime T <C |3Ef ,-,„p(0)|, we obtain 
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where Aq is the spectroscopic gap, and /i = 
{\ / ^())\dA{ip) I dip\^^^^^ is the opening rate of the gap 
function at the node. In analogy with the thermal con- 
ductivity [41, the remaining effect alT ^T* of inelastic 
scattering is a modification of the T ^ Q asymptotic. 



When we write explicitly J^f 



and ^ = v^/(l + Am) in Eq. (|3}, one factor 1 + Air, 
remains in the denominator. Within the bare theory, this 
result can be traced back to that the spectroscopic gap 
Ao in the weak-coupling limit is replaced by the strong 
coupling off-diagonal function W(0). 

In summary, we have discussed the temperature- 
dependence of the thermoelectric effect that results from 
the interplay of inelastic electron-boson scattering and 
elastic impurity scattering. At high temperatures, T > 
T*, particle-hole symmetric inelastic scattering domi- 
nates and the thermoelectric coefficient is suppressed. On 
the other hand, at low temperatures, T < T*, particle- 
hole asymmetric elastic impurity scattering dominates 
and the thermoelectric coefficient is enhanced. 
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